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Abstract. We study the propagation of waves from the photosphere to the chromosphere of 
sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and 
polarimetric spectra of Si I A 1082.7 nm and He I A 1083.0 nm we retrieve the line-of-sight velocity 
at several heights. The analysis of the phase difference and amplification spectra shows standing 
waves for frequencies below 4 mHz and propagating waves for higher frequencies, and allows us 
to infer the temperature and height where the lines are formed. Using these observational data, 
we have constructed a model of sunspot, and we have introduced the velocity measured with the 
photospheric Si I A 1082.7 nm line as a driver. The numerically propagated wave pattern fits 
reasonably well with the observed using the lines formed at higher layers, and the simulations 
reproduce many of the observed features. The observed waves are slow MHD waves propagating 
longitudinally along field lines. 



1. Introduction 

The magnetic field of the sunspots affects the propagation of waves and produces the appearance 
of new wave modes not present in the quiet Sun. The study of waves in sunspots is interesting 
from several points of view. It gives information about the atmospheric structure and dynamics. 
Also, waves are believed to constitute one of the heating mechanisms of the chromosphere. 
Observations show that the properties of waves greatly vary with the region where they are 
observed: at the umbral photosphere they present 5 minute period with amplitudes around 
100 m s _1 , while the chromospheric umbral oscillations have basically 3 minute period and 
amplitudes of a few kilometers per second. 

Recent simulations have been able to reproduce many of the observed properties of waves 
between the photosphere and the chromosphere by means of the introduction of a photospheric 
pulse which drives a wave spectrum close to the solar one [lj . In this work we go a step further and 
we have developed numerical simulations where the real photospheric velocity pattern retrieved 
from observations in the Si I A 1082.7 nm line is used to drive the simulated wave pattern. 
The oscillations obtained with the simulation at the high photosphere and chromosphere are 
compared with those measured with lines which are formed at those heights. 



2. Observations 

We have obtained simultaneous cospatial maps of the LOS velocity measured with photospheric 
lines (Si I A 1082.7 nm and the Fe I lines at 396.54, 396.61, 396.66, 396.74 and 396.93 nm from 
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Figure 1. Spectra of the Call H line. The 
lines used and their formation height are 
marked. 
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Figure 2. Spectra of the region at 1083 nm. 
The lines used and their formation height are 
marked. 



the wings of Call H) and chromospheric lines (Ca II H core and He I A 1083.0 nm). The phase 
and amplification spectra between several pairs of lines in this set of observations have been 
analized before in [2j. From the fit of the observations to a simple model of linear upward wave 
propagation we have infered temperature and formation height of the lines. Figures [1] and El 
show the observed spectral region. The relative formation height of the spectral lines is written 
next to each line. 

3. Numerical simulations 

The numerical simulations used in this work were performed with a 3D non-linear MHD code [1] . 
Energy losses are implemented following Newton's cooling law, accounting for the damping of 
the temperature fluctuations due to radiative losses with Spiegel formula for the typical cooling 
time [3]. 

We have constructed a MHS model of sunspot using observational data. The stratification 
of the thermodynamical magnitudes, obtained from the inversion of the umbra and the quiet 
Sun photospheric spectra of Si I line (using SIR [I]) are set at the axis of the sunspot and 
at 40 Mm far from it, respectively. These two atmospheres are smoothly merged along the 
radial distance which separates them [5]. The model is an axysymmetric thick flux tube with 
distributed currents and has no twist. In the horizontal directions the computational domain 
covers 8.4 x 5.2 Mm, with an horizontal spatial steps of Ax = Ay = 100 km. In the vertical 
direction it spans form z = —0.6 Mm to z = 1 Mm, excluding the special layer added at the 
boundaries to absorb the waves without reflections, with a spatial step of Az = 25 km. The 
MHS model of the sunspot is pertubed around the photospheric formation height of the Si I line 
with the velocity measured with this line. This driver generates an oscillatory pattern identical 
to the observed one at the photosphere, and these waves propagate upward to higher layers. 

4. Results 

In order to compare the numerical simulation with the observational data, we have assigned a 
fixed z to the formation height of each spectral line, and we assume that the vertical velocity at 
that height corresponds to the velocity measured from the Doppler shift of the line, (see Figures 
□ and[2j) [21 . 

Figure [3] shows the observational and numerical velocity maps, corresponding to the 
chromosphere. It is the highest spectral line observed, and these waves have propagated upward 
about 800 km from the formation height of the Si I line in order to reach this layer. Note 
that in the simulated velocity map (middle panel) the velocity signal is almost nule during the 
first 2 minutes, due to the time spent by the slow waves to cover the distance between the 
driver and this height travelling at the sound speed. During this travel the period of the waves 
is reduced to around 3 minutes and their amplitude increases, reaching peak-to-peak values of 
almost 8 km s _1 . Bottom panel of Figure [3] shows that the oscillations develop into shocks. This 
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Figure 3. Velocity maps for 
the chromosphere inside the umbra. 
Top: observational (measured from 
the Doppler shift of the He I line); 
middle: numerical (vertical velocity 
at the formation height of the He i) ; 
bottom: temporal variations of 
numerical (blue) and observed (red) 
velocities at x = 0.9 Mm. 
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Figure 4. Power spectra of the 
velocity measured with the He I line 
(red dashed line) and the simulated 
at the corresponding height (blue 
line) . 



behavior is well reproduced by the numerical simulation. The simulated velocity map reproduces 
reasonably well the observed oscillatory pattern. Only in the temporal lapse between t = 27 
and t — 40 min the simulated pattern differs signiflcatively from the observations. 

The power spectra at chromospheric heights is dominated by 3 minute oscillations in the band 
between 5 and 8 mHz (Figure HJ. The observational power has three power peaks in this band, 
located at 5.5, 6 and 7 mHz. The simulated power is concentrated at a single peak between 
the two highest peaks of the observations. The simulated peak at 5.5 mHz is lower than the 
observed one. The simulations also reproduce the power peaks at 7.7 mHz and 9 mHz, and the 
low power at frequencies below 5 mHz. At frequencies above 13 mHz the simulated power is 
slightly higher than the observational one. 

Phase spectra give the delay between two signals at every frequency. The phase difference 
between the velocity signals measured in the photospheric Si I line and the chromospheric He I line 
is zero for frequencies below 4 mHz (left panel of Figure [5|). At these frequencies the atmosphere 
oscillates as a whole, i.e., the waves are stationary. From 4 mHz to 7 mHz, the phase difference 
increases linearly with the frequency. It indicates that waves at these frequencies propagate from 
the photospheric layer, where the Si I line forms, to the chromospheric layer, where He I forms. 
Note the excellent agreement between observational and numerical data from to 7 mHz, where 
the coherence of the observations is high. At higher frequencies the coherence of the observed 
phase difference is lower, but the simulated one keeps its linear increase. With regards to the 
amplification spectra, for frequencies above 1.5 mHz the simulated spectra reproduces properly 
the observed one. 
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Figure 5. Left panel: Phase difference (top) and amplification spectra (bottom) between Si I 
and He I velocities. Right panel: Phase difference (top) and amplification spectra (bottom) 
between Si I and Fel velocities. The background color in the top panels indicates histograms of 
the relative occurrence of a given value in the observed phase (from black -low- to yellow -high-). 
Green crosses are the results of the simulation for all the spatial points. The white line is the fit 
to a model of linear wave propagation [2j. The bottom panels show the ratio of the power at the 
heights of the two signals (amplification spectra) as a function of frequency for the observation 
(black solid line), the simulation (black dashed line) and the model (red line). 



Right panel of Figure [5] shows the phase and amplification spectra between two photospheric 
lines, the Si I and the Fel A 396.93 lines. The observed phase difference has high coherence 
between 2 and 8 mHz, and in this frequency range the simulated phase delay fits the observational 
one reasonably well, showing a Acj) = for frequencies below 4 mHz and a slow increase for higher 
frequencies. For frequencies below 2 mHz and above 8 mHz the observed phase difference spreads 
out and has lower coherence. The behaviour of the simulated amplification spectra is similar to 
the observed one between and 8 mHz, but the numerical amplification is significatively higher 
at the peak around 6.5 mHz. 

5. Conclusions 

We have studied the propagation of waves in sunspots from the photosphere to the chromosphere, 
using observationally driven numerical simulations. We have been able to simulate the travel 
of the waves through the sunspot atmosphere, with a remarkable match between observations 
and simulations. The numerical calculations reproduce the wave pattern (Figure the power 
spectra (Figured]) and the phase and amplification spectra between several pairs of lines (Figure 
[5]), indicating that observed waves are slow longitudinal acoustic waves dampled by radiative 
losses. 
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